Magnetic field effects on spin relaxation in heterostructures 
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Effect of magnetic field on electron spin relaxation in quantum wells is studied theoretically. 
We have shown that Larmor effect and cyclotron motion of carriers can either jointly suppress 
D'yakonov-Perel' spin relaxation or compensate each other. The spin relaxation rates tensor is 
derived for any given direction of the external field and arbitrary ratio of bulk and structural 
contributions to spin splitting. Our results are applied to the experiments on electron spin resonance 
in SiGe heterostructures, and enable us to extract spin splitting value for such quantum wells. 



I. INTRODUCTION 

Spin dynamics of charge carriers in semiconductors 
and semiconductor nanostructures attract much atten- 
tion during the last decade. Theoretical and experimen- 
tal efforts are made in field of creation, conservation and 
manipulation of spin polarized electrons. In this regard 
the control of spin dephasing processes is a task of prime 
importance. 

The main spin relaxation mechanism in high qual- 
ity semiconductor quantum wells is that proposed by 
D'yakonov and Perel'i^: electrons lose their initial spin 
during free motion between scattering events due to 
wavevector dependent effective magnetic field i^ This 
field is determined by the conduction band spin split- 
ting in the systems lacking an inversion center. In bulk 
zinc-blende lattice semiconductors the spin splitting is 
cubic in electron wavevector k (Dresselhaus term), and 
in quantum wells it can be presented as a sum of linear 
and cubic contributionsi^ Possible asymmetry of a quan- 
tum well gives rise to an additional linear in wavevector 
contribution to the spin splitting (Rashba term)i^ These 
two bulk and structural asymmetry terms have different 
symmetrical properties and lead to the giant spin relax- 
ation anisotropy in the heterostructure planed. The simi- 
lar situation is observed in SiGe quantum wells: although 
bulk Si and Ge possess inversion center, quantum wells 
grown from these materials can lack such a center and 
allow the spin splitting of the electronic subbands which 
can have symmetry of Dresselhaus or Rashba terms. ^ 

One of the main features of D'yakonov-Perel' spin 
relaxation mechanism is that it can be suppressed by 
all processes which change electron wavevector such 
as scattering from impurities, interface microroughness, 
phonons, electron-electron collisions^ and because of cy- 
clotron motion of carriers in the external magnetic fieldi^ 
Besides, magnetic field slows down spin relaxation due 
to Larmor spin precession around the external field 
directioniSiS In Ref. [l^ it was shown that depending on 
type of linear in k terms, the interference between the 
Larmor and cyclotron effects will result in either sup- 
pression of spin relaxation due to joint contributions of 
both effects or to their mutual compensation leading to 
partial or full recovery of spin relaxation rate. However, 
expressions presented in Ref. llfl were obtained in the 
limit when only one contribution to the spin splitting is 



present, moreover as we will see below they predict in- 
correct dependence of the transverse relaxation time T2 
on magnetic field orientation. 

The aim of present paper is to develop a consistent the- 
ory of D'yakonov-Perel' spin relaxation in classical mag- 
netic fields in presence of both Dresselhaus and Rashba 
splittings. In contrast to Ref. [13 the expressions for all 
components of spin relaxation tensor will be obtained. 
When both Dresselhaus and Rashba terms are compa- 
rable in magnitude the spin relaxation in the quantum 
well plane is strongly anisotropic and correct definitions 
of longitudinal Ti and transverse T2 spin relaxation times 
is required. We present an expression for spin susceptibil- 
ity with allowance for the relaxation anisotropy and show 
how electron spin resonance (ESR) linewidth is connected 
with spin relaxation tensor components. We analyse re- 
cent experiments on ESR in SiGe quantum wells^" and 
extract electron spin splitting from the angular depen- 
dence of the linewidth. 



II. THEORY 

Spin-dependent contribution Tiso hr the effective 
Hamiltonian consists of two components. First, in 
the asymmetrical quantum wells there exists so-called 
Rashba term^ 
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where constant a is nonzero due to the heteropotential 
asymmetry, Ui {i = x,y,z) are Pauli matrices. Axes x 
and y are directed along [100] and [010] respectively. 
Second contribution has a Dresselhaus formi 
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Piio-yky - a^kx) + f3i(<7xkxkl - aykykl). (2) 



In the quantum wells grown from zinc-blende lattice ma- 
terials P3 is the bulk spin-orbit constant i2i4 Even if 
(3^ = (like in SiGe quantum wells) /3i is nonzero due to 
the asymmetry of chemical bonds at heterointerfaces&ii 
Sum of Rashba (Q) and Dresselhaus-like Q terms can 
be conveniently represented as 



T^so — 'Hd + Ti-R 



n. 
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where components of vector a are Pauli matrices and 
rjfc is the effective Larmor frequency of spin precession 



in the spin-splitting induced field. One can see that vec- 
tor rjfc contains both first and third harmonics of the 
angle between the wavevector k and x-axis. Spin preces- 
sion around ilfc results in spin relaxation while processes 
changing electron wavevector such as scattering and cy- 
clotron motion slow it down. 

The application of the external magnetic field B leads 
to (a) cyclotron motion of carriers in the quantum well 
plane, its frequency is given by 
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where B^ is normal component of magnetic field, e and c 
are elementary charge and light velocity, m is the electron 
effective mass and (b) Larmor precession of electron spins 
with angular velocity u:i^. The components of Wl can be 
presented as 



WL, 



^ R 



(5) 



where ^b = eh/2moc is Bohr magneton, ttiq is free elec- 
tron mass and gij are the components of the ^-factor 
tensor for electrons in the quantum well, and summation 
over the repeated indices is omitted. The direction of the 
vector wl defines the direction of Larmor precession axis. 
In what follows it is convenient to introduce vector wc 
directed along the growth axis with absolute value being 
equal to the cyclotron frequency. 

We will use density matrix formalism in order to de- 
scribe spin relaxation. Electron spin density matrix can 
be decomposed as 

Pk = fk + Sk ■ a- , 

where fk is the spin-averaged occupation of the fc state 
and Sk is the spin vector for electrons in this state. Ki- 
netic equation for spin distribution Sk can be written in 
the following formL^ 



QSl, 

— — + Sfc X (wl + ^k) 



Ask + Q{sk}^0, (6) 



where the second term describes spin precession in the 
presence of both external field and spin splitting in- 
duced effective magnetic field Eq. (PJ, the third term 
comes from wavevector variations under cyclotron mo- 
tion with operator A defined according to {Ask)i = 
u>c[fc X dsk^i/dk], and the last term is the collision inte- 
gral. For simplicity we will consider only elastic scatter- 
ing processes and describe them using two times ri and 
T3, responsible for relaxation of first and third harmon- 
ics of distribution function;- Eq. ^ is valid for classical 
magnetic fields, quantum effects on spin relaxation were 
discussed in Refs. Iiall3l for strong and weak fields re- 
spectively. 

Considering spin-splitting as a small perturbation 
(fifcTi, QkTs <C 1) one can present spin distribution func- 
tion Sfc as 



where s^ is quasi-equilibrium distribution function and 
Ssk is a non-equilibrium correction arising due to spin 
splitting. After summation © over the wavevectors we 
arrive to the balance equation describing relaxation of 
the total spin S = J2k ^k 



dS 
'dt 



5 X u;l + rS' = 0, 



(7) 



where the tensor of inverse spin relaxation times T 
defined according to 
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6sk X Ofc, 



(8) 



and the non-equilibrium correction 5sk should be deter- 
mined from the solution of Eq. I^. 

We will be interested in the case of strong enough mag- 
netic field : f^lri, flj^r^ <C wl- Then all non-diagonal 
components of the tensor T associating spin relaxation 
along external field and in the transverse plane are zero. 
Noting that the first and third harmonics of spin distri- 
bution have independent dynamics we present T in the 
form 



f = fW+f(3), 



(9) 



where the upper index denotes the number of angular 
harmonic. 
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FIG. 1: Frames of axes, z \\ [001] is a growth axis. Axes 
xi 11 [110] and yi \\ [110]. The frame of axis x', y' and z' is 
connected with the external magnetic field: z' \\ ujl, y' is a 
perpendicular to o^l lying in a quantum well plane, and x' is 
chosen so that x' , y' and z' form a right-hand triple. 



It is convenient to present components of spin relax- 
ation tensor in the new frame of coordinates with z' axis 
directed along Larmor frequency vector wl, and y' axis 
lying in the quantum well plane and being perpendicular 
to the ujl (Fig. P). Let 9 be the angle between z' and z 
and ip be the angle between y' and yi \\ [110]. One can 
show that for electrons with fixed energy E the non-zero 
components of tensors F^^^ F'^^ can be written as 
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Here k = y/2mE/h, i — x',y', /3 = (3i — P^k'^ / A is a. coefficient at ffist harmonics in Dresselhaus splitting ^, 
Dj^± = 1 + (ojl ± J^c)^Tj, where j = 1, 3, and 



Ni^^>^> = a^ [l + cos^ e + {ujc+ WL cos 0) V^ + {ujc cos 6* + wl) ^i ] 
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+/3^ [l + cos2 6* + (cjc - t^L COS 6')^t2 + (wc cos 6* - wl)^t2] + 2a/3 cos 2^; sin2 e[l + (cj^ + '^l)'^? ] : 

iVi^x'x' = («'+/?') {1 + sin2 0(l +cc;2r2)[l + {lj^ + iul)T^)]{Di^+Di^_)-^} + 2a/3cos2V' , 

Ni^y'y' = (a^ +/3^) {1 - sin2 6'[l + {luI - 3luI)t^]{Di^+Di,^)-^ujItI} ~ 2a/3cos2V' , 



N^^z'z' = 1 - sin" - <! 1 



cos^ 1 + (o^L + 9cjq)t3 1 — 6uJi,uJcTi 
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Eqs. pi)|) and (|ll|l describe spin relaxation of two-dimensional electron gas subject to classical magnetic field. These 
formulae are valid for arbitrary ratio of Dresselhaus and Rashba contributions and for arbitrary orientation of the 
external field. Generalization of Eqs. (|10|l . (|ll|l for the case of any given energy distribution of electrons can be 
carried out in the standard wayi 



III. RESULTS AND DISCUSSION 



A. Longitudinal relaxation 



We will start the analysis of the expressions Eqs. (|1(J|I and pi|) for the spin relaxation tensor components in the 
case when magnetic field is normal to the quantum well structure, i.e. 6 = Q. Then the expression for the longitudinal 



spin relaxation rate 1/Ti — T ^ 
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Eq. H12(l shows that external magnetic field has differ- 
ent effects on spin relaxation depending on the dominant 
type of spin splitting. When Rashba term is dominant 
(/33 = 0, /3 = 0, a ^ 0) and electron 5-factor is positive 
the Larmor and cyclotron effects partially compensate 
each other. In the case when first harmonics in Dres- 
selhaus splitting is dominant (/3 S> ii-ik'^ ■, a = 0) both 
Larmor and cyclotron effects jointly slow down spin re- 



laxation. If the spin splitting is dominated by the third 
harmonics then Larmor and cyclotron effects compensate 
each other. 

Such interference behavior is a result of the fact that 
according to Eq. ((SJ spin relaxation rate is determined 
by the vector product of nonequilibrium correction (Js/, 
and effective spin precession frequency $7^. If no mag- 
netic field is applied then bsy, is orthogonal to $7^ and 



03 



1.2 



{ 1.0 















"~~~- 


^^ 




^' 




^\ 








\. 


" t 


~~~~-_ ''^ 





--''' 


e(deg) 


■ 





30 60 

e(deg) 



90 



3 




30 60 

e(deg) 



90 



FIG. 2: Longitudinal spin relaxation rate vs. magnetic field 
tilt angle 6 for SiGe based quantum well. Spin relaxation rate 
is measured in the units 2T\k^5^ /h^ where 5 is a coefficient at 
the first harmonics of spin splitting. Solid curve corresponds 
to Rashba term (a = 5, /? = 0), dashed curve to Rashba term 
{a = Q, 13 = 5). The contribution of the third harmonics is 
neglected. At ^ = the parameter loqti — 1. The parameters 
of calculation were taken from Ref. Jj; g — 2, m — 0.196mo. 
Inset shows nonmonotonous behavior of 1/Ti{9) when Dressl- 
haus term dominates. 



spin relaxation rate is fastest. External magnetic field 
results in the precession of in-plane spin component Ss^ 
around Larmor vector wl and in cyclotron rotation of the 
rife. If these rotations are not synchronous then vectors 
6sk and ilfc are no longer orthogonal and spin relaxation 
rate becomes slower. In the case of synchronous rotation 
of in-plane spin and effective magnetic field the spin re- 
laxation is not suppressed. It is important to note that 
the direction of cyclotron rotation of $7^ depends on the 
dominant spin splitting mechanism. If Rashba term is 
largest then effective magnetic field Jlfc is perpendicular 
to the wavevector and rotates in the same direction as 
k. Thus, spin relaxation becomes slower by the factor of 
1 + (cjl — (^c)'^Ti and Larmor and cyclotron effects par- 
tially compensate each other provided electron g-factor 
being positive. If spin splitting is due to the first harmon- 
ics of Dresselhaus term then cyclotron rotation of ilfc and 
Larmor precession at positive (7-factor have opposite di- 
rections and slow down spin relaxation by the factor of 
1 -I- (wl + (^c)'^Ti. When third harmonics dominates, the 
rotation of fife and Larmor precession take place in the 
same direction but the rotation frequency of flk is three 
times higher than cyclotron one. 

Figure|5]shows longitudinal spin relaxation rate 1/Ti = 
Tz'z' dependence on the magnetic field tilt angle 6 calcu- 
lated for SiGe structure. The case of dominant Rashba 
term is shown by a solid line, and dashed line presents 
the results when spin splitting is determined by Dres- 
selhaus term. The contribution of the third harmonics 
is neglected in all calculations. The qualitatively differ- 
ent angular dependencies of 1/Ti are clearly seen. When 



FIG. 3: Longitudinal spin relaxation rate vs. magnetic field 
tilt angle in the case of equal Dresselhaus and Rashba 
terms(Q: — /3) for different orientations of magnetic field. Spin 
relaxation rate is presented in units 2T\k^o? /h^ . The values of 
parameters used in the calculation are the same as for Fig. |21 



Rashba contribution is a main one the spin relaxation 
rate decreases with increase of in-plane component of 
magnetic field because tilting of the field results in the 
decrease of the cyclotron frequency wc- When Dressel- 
haus term is dominant the dependence \/Ti{9) is non- 
monotonous (see inset): first, spin relaxation slows down 
owing to the decrease of projection of Ofc on the plane 
perpendicular to the magnetic field. At larger 9 spin re- 
laxation accelerates due to cyclotron frequency decrease. 

Figure |31 shows longitudinal spin relaxation rate for 
different orientations of magnetic field in the case when 
Rashba and Dresselhaus terms have equal strengths (a = 
(3). Different curves correspond to different angles be- 
tween magnetic field and main axes of the structure. One 
can see that spin relaxation rate is highest aX ip = Q, i.e. 
for magnetic field inclined to the xi \\ [110] axis and is 
lowest sX ij} — -K 12 (magnetic field inclined to y\ || [110]). 
This result is a direct consequence of the interference of 
linear in fc spin dependent terms resulting in the sup- 
pression at a = /3 of the spin splitting along [lIO].^ For 
the same reason spin relaxation rate calculated with al- 
lowance for the first harmonics of the wavevector only 
goes to zero when magnetic field is directed along [110], 
the account for the third harmonics in spin splitting leads 
to the finite spin lifetime along j/i^ 

Different behavior of spin relaxation rate for different 
values of ■0 can be explained taking into account that 
in normal magnetic field spin relaxation is suppressed 
mostly by cyclotron effect since oji^/loq — gm/2mQ « 0.2. 
With increase of 9 cyclotron frequency decreases, and if 
— 7r/4 < ip < IT /A effective spin splitting increases. There- 
fore spin relaxation becomes faster if — 7r/4 < tp < 7r/4. 
At -0 = 7r/4 and ojcti — 1 spin relaxation rate should not 
depend of tilt angle if one completely neglects Larmor ef- 
fect. With tilting of magnetic field cyclotron frequency 
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FIG. 4: Transverse relaxation time vs. magnetic field tilt angle 6 for SiGe structure, (a) Open squares show ESR linewidth 
measured in Ref. HoL Solid curve presents the results of calculation when Rashba term dominates. Parameters loc ~ 2.96 x 
10^^ s~^ at normal magnetic field orientation, ojl = 5.93 x 10^" s~^, and momentum relaxation time ri = 10~^^ s were taken 
from experiment. Spin precession frequency is Q.kp = akp/h = 4 x 10* s~^. (b) Dresselhaus and Rashba terms are equal 
{a — 13). Spin relaxation rate is measured in the units of nfc a jfi . Parameters of calculation are the same as for Fig. |2K- 



decreases and Larmor effect becomes more important re- 
sulting in slight decrease of relevant curve (V' = 7r/4). 
If 7r/4 < ip < 37r/4 then competition between decrease 
of wc, decrease of effective spin splitting and Larmor ef- 
fect is possible leading to nonmonotonous angular depen- 
dence of the spin relaxation rate. However at loqTi — 1 
taken for calculation presented in Figl^J) spin relaxation 
slows down. 



B. Transverse relaxation 

Now we turn to the analysis of spin relaxation in the 
plane perpendicular to the external field. In the ESR 
experiments the spin susceptibility is measured which in 
the limit of low microwave power has a following form 
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Xoc 



WL 



i/T2 



(13) 



where to is microwave frequency and T2 is a transverse 
relaxation time. Eq. \V6\ is valid when microwave fre- 
quency is close to Larmor one, small renormalization of 
resonance frequency ^ T2 is neglected. Transverse re- 
laxation time can be expressed through the components 
of spin relaxation tensor as 



1/^2 = 



r 



y'y' 



(14) 



and ESR half-maximum full width Auj = 2/r2. Eq. ^ 
generalizes results of Refs. 10 14 for the case of arbitrary 
orientation of magnetic field and anisotropic spin relax- 
ation. One can see that the combination T^'x' + '^y'y' is 
invariant to the choice of axes in the transverse plane, 
therefore resonant peak width is independent of mi- 
crowave field polarization. 



If the external field is oriented along the quantum well 
growth axis then I/T2 has the form 



l/T2 = 
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In the particular case of /3 = and loq — this expression 
coincides with equations obtained in Ref. Il4l However 
it is in sharp disagreement with results of Ref. llCt where 
zero spin relaxation rate is predicted for normal mag- 
netic field. Those results were based on the assumption 
that transverse spin relaxation is governed by the normal 
component of J7fc. This assumption is invalid, since spin 
relaxation along, say, x' is governed by both ^k,y' and 

Figure ^ shows the results of calculation of the trans- 
verse spin relaxation rate I/T2 in SiGe structure versus 
external magnetic field orientation. The transverse re- 
laxation rate monotonously increase with tilt angle (solid 
curve on Fig. 0^). It comes from the fact that spin relax- 
ation in the plane normal to magnetic field is suppressed 
owing to cyclotron effect, see Eqs. (tTUI) and (fTC|l . Cy- 
clotron frequency is maximal if magnetic field is normal 
to the quantum well plane and decreases proportionally 
to cos 9 with tilting magnetic field. 

Open squares on Fig. ^ show the ESR linewidth mea- 
sured in SiGe quantum well structurei^S The parame- 
ters of calculation were taken equal to the experimental 
ones: ujc = 2.96 x 10^^ s~^ at normal magnetic field, 
cjl — 5.93 X 10^° s~^, the momentum relaxation time 
Ti — 10^^^ s. We used the effective frequency of spin 
precession on the Fermi level ^kp = akp/h as a fitting 
parameter. The best agreement is at ^kp = 4 x 10* s~^. 
The order of magnitude of flkp agrees with estimations 
of Ref. [Ij. Obtained splittin g is lower than value 
6.3 x 10® s""'^ presented in Ref. (I^, because fitting in 



that paper was done with incorrect dependence \/T2{0). 

In the case when Dresselhaus and Rashba terms are 
equal in magnitude (a = /3), the dependence of trans- 
verse relaxation rate 1/T2{0) on magnetic field tilt angle 
is qualitatively the same. The slope depends on the in- 
plane rotation angle ^p reflecting anisotropy of the sys- 
tem. 

In conclusion, we have studied in detail effects of mag- 
netic field on electron spin relaxation in SiGe quantum 
wells. We have shown that, depending on electron g- 
factor sign and dominant spin splitting mechanism, Lar- 
nior and cyclotron effects are either compete with each 
other or jointly suppress spin relaxation. The expression 



for ESR linewidth was derived for the case of anisotropic 
transverse relaxation. Fitting by the derived expressions 
the experimental data we obtained spin splitting value 
for SiGe structure studied in Ref. llOl 
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